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Abstract. On the bases of the improved approximation for the spectral function of one- 
particle states the Landau-Silin kinetic equations for the normal Fermi liquids of neutral 
and electrically charged particles are shown to be valid at finite temperature above the 
temperature of superfluid transition. 

1. Introduction 

The exclusive successfulness of the phenomenological Landau normal Fermi liquid theory 
[1] in predicting and describing a set of new phenomena, among them the zero sound in liquid 
He 3 and spin waves in nonferromagnetic metals, made this theory a subject of investigation 
on the basis of strict microscopic theory. Most attention was devoted to the derivation of 
the kinetic equation for the quasiparticle distribution function. The initial derivation of 
this equation was produced by Kadanoff and Byam [2] on the basis of the quasiparticle 
(QP) approximation for the spectral function and was continued by some followers who 
used so called extended quasiparticle (EQP) approximation [S]-[7j. However, in all these 
cases the second Poisson bracket in the left side of Kadanoff-Baym (KB) generalized kinetic 
equation could not be eliminated in a lawful mathematical way (see below). This fact made 
the temperature range of validity of the kinetic equation for quasiparticle distribution very 
narrow, strictly speaking the theory was proved to be valid only in the vicinity of absolute 
zero. 

Experimental discovery of the superfluidity of He 3 at the temperature lower than the 
temperature at which the zero-sound in the normal Fermi liquid was discovered, and theo- 
retical works devoted to the description of the superfluid state [HI E] , left no room for the 
temperature range of validity of the equation for normal Fermi liquid. At the same time 
theoretical predictions of this theory turned to be in a perfect numerical agreement with 
the experimental data. This fact stimulated a second wave of attempts to deriving the ki- 
netic equation for normal Fermi liquid, but as it was mentioned above, the result was not 
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achieved to a satisfactory extent. Further development of KB theory went in the direction 
of the constructing of a nonlocal quasiparticle kinetic equations [TUl [TT] . a development of 
a detailed selfconsistent microscopic treatment of arbitrary initial correlations in the sys- 
tem [IS], etc., but the question of the temperature range of validity of the Landau's kinetic 
equation remained open. 

2. Spectral function and kinetic equation for normal Fermi Liquid 

The problem consisted in a mathematically lawful elimination of the second generalized 
Poisson bracket (a" puzzling term" in the terminology accepted in [5]) in the generalized 
Kadanoff-Byam kinetic equation written in the collisionless approximation in the case of 
slowly varying in space and time disturbances when only the first derivatives with respect 
to T and R are taken into account [2]: 

(1) [u - e(pu; RT),g < (pcu; RT)} + [Re g(pu; RT), a < (pto; RT)} = 0, 



here [A,B]-the generalized Poisson bracket, defined by the expression: 

dA dB dA dB 
dudf " dfdu 



(2) [A B] = —— - —— - V,A . V^B + V^A . VpB, 



e(pu; RT)-the energy of a particle, defined by the equality: 

(3) e(pcu; RT) = E HF (p; RT) + Re a(pu; RT), 

E HF is a one-particle energy in the Hartree-Fock approximation, Re a is a real part of the 
correlation energy of a particle related to the imaginary part through the Hilbert transform, 

n M - , „ f°° duj'T{puj';RT) 

(4) Re a(puj- RT) = P KF / . 

Here, P refers to a principle value integration. 

It was shown in [2] that the function g entering the Eq. (1) can be taken in the form 

(5) g(pZ; RT) = [Z- E HF (p; RT) - Re a(pZ; RT)}' 1 . 

Eqs. (4) and (5) lead to the following general form for a spectral function of one-particle 
states in the system: 

T(pu;RT) 



(6) a(pu; RT) 



[u - e(pu;RT)] 2 + 



V 2 (pw;RT) 
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The spectral function satisfies the exact sum rule: 

(7) f ^a(pu; RT) = 1. 

Correlation function g < is related to a spectral function by the equality [2]: 

(8) g<(puj; RT) = a{puo; RT)f(puj; RT). 

The complexity of the expression (6) makes to search for a certain approximations for the 
spectral function which could be successfully used in calculations. 

Quasiparticle approximation for the spectral function (6) corresponds to the case T — > 0, 
when the quasiparticles are stable: 

(9) a QP = 2nZ(p; RT)5[u - E(p; RT)}, 
where E{p; RT) is the solution of the equation: 

(10) E(p; RT) = E HF (p; RT) + Re a(pu; RT)l =mRr) , 
and the renormalizing factor Z(p; RT) is given by the expression: 

(11) Z-\p-RT) = \ 



... , Fyr ,, i d Re a(ptu; RT) 



® U oj=E{p;RT) 

Using the approximation (9) one follows the KB way to deriving the Landau's kinetic equa- 
tion when the second Poisson bracket in Eq. (1) is dropped considering o r< = fT to be 
negligibly small what is valid only at u = \i at zero temperature. 

The extended quasiparticle (EQP) [3]-[7j and improved extended quasiparticle (iEQP) [1] 
approximations for the spectral function are written in the form: 

r 

(12) a EQP = 2nZ5{u - E) + P- 



uj-Ef' 

r 



(13) a lEQP = 2nZ5(uj - E) + ZP- 

[oj — hiY 

where E and Z are determined by the Eqs. (10) and (11) correspondingly. 

It is easy to verify that the approximations (12) and (13) do not satisfy the KB equation 
for a nonequilibrium spectral function: 

(14) [u - e{poo; RT), a(pu; RT)} + [Re g(pu; RT),T(pu; RT)} = 0, 
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and do not lead to the elimination of the second Poisson bracket in Eq. (1) in the case of 
finite T. Only such a mathematically lawful elimination of this term for finite values of T 
would testify the validity of the Landau's equation at a finite range of temperature. 

The expressions (12) and (13) both were obtained on the basis of the general expression 
(6) by means of the Taylor expansion in powers of V in the frame of different approximations 
[3]-[7]. As a result, some important factors were lost. In reality there does not exist a 
mathematically strict correct form for the expansion of (6) in power series of T which starts 
with the delta function when T = 0. 

Another expression for the spectral function can be offered on the basis of the following 
consideration. We utilize a well known relation of the Fourier transform in the case of a 
constant value of c: 

(15) / e-^ c e ltx dt = - , c> 0. 

J-oo C + X 2 

Expanding the first exponent in the left side of (15) in Taylor series, we get: 

2r f°° r 2 t 2 

(16) -^— 2 = I (l-c\t\ + C -^-...)e^dt 



c 2 + x 2 1 1 2! 



Now we use the formulas equivalent to those represented in [13,14]: 

POO 

(17) / t {2n) e ttx dt = 2n(-t) 2n 5 {2n \x), n = 0,1,2, 



J — oo 
oo 



;i8) / \tf n+1 h ltx dt = -2sin((2n + l)J)(2n + l)! ' 



2 |x| 2n+2 ' 



n = 0,1,2, 



If the quantity V in Eq. (6) was not a function of u, then the expressions (17), (18) 
would lead to a strict correct expansion of the spectral function (6) in terms of the power 
series of T. In the case of T depending on oo, one can rely only on the first two terms of 
the expansion: the delta function independent of T and the term proportional to T. Then, 
taking into account that 

(19) uj - e(pu; RT) = Z' 1 ^; RT)[uj - E(p; RT)], 
we come to the following approximation for the spectral function (6): 

r 

(20) a = 2-nZ5(u -E) + Z 2 P 



(lu-E) 2 ' 
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It is not difficult to verify that the approximation (20) satisfies the sum rule (7), with the 
same precision as approximations (12) and (13) do, but contrary to them the approximation 

(20) satisfies the Eq. (14) for the spectral function and eliminates the second Poisson bracket 
in the left side of Eq. (1). Indeed, when we substitute Eq. (20) into Eq. (14), the first term 
in the right of this expression gives: 

(21) [u - e, 2nZ5(u - E)} = 2tvZ[(uj - E)Z~ 1 , 5(u - E)} = 0, 
due to the property of the generalized Poisson bracket 



[AJ(A)]=0. 



The second term in the right side of Eq. (20) leads to the expression: 



(22) 



to — e, Z z 



1 



(w - E)' 



:[U-E,T], 



{u-E) 2 . 

the second poisson bracket in Eq. (14) due to Eq. (5) gives the expression: 

1 



(23) 



[Reg,r] = Z 



l(u-E) 



r 



Thus, the ansatz (20) satisfies the Eq. (14) exactly. Now we consider the Eq. (1). Using 
Eq. (8), we get: 



(24) 



[u-e,af] + [Reg, fT] = 0. 



The substitution of the first term in Eq. (20) into Eq. (24) leads directly to the kinetic 
equation for the quasiparticle distribution function n(p; RT) [2]: 

dn 



(25) 



0T +VpE. V A n-V A E. V p -n = 0, 



n(p;RT) = f(puj;RT) \ u=mAr) ■ 
The second term in Eq. (20), being substituted to Eq. (24), gives the expression: 



(26) 



u — e, 



(uj-E) 2 \ 



Z[uj-E,Tf] 



(uj-Ef 
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which is eliminated by the second poisson bracket in Eq. (24): 

[Reg,a < } = \^-^,rf} = Z 
(27) z^-E]-!—. 

Thus, the kinetic Eq. (25) is valid for finite values of T and correspondingly, for nonzero 
temperatures. For qualitative estimation of the precision, we can use the third term in the 
expansion (16) which is proportional to T 2 . Substituting this term into Eq. (24), it is not 
difficult to show that this equation is valid up to the terms of order T 2 . 



uj-E 



r/ 



3. KINETIC EQUATION FOR NORMAL FERMI LIQUID IN A MAGNETIC FIELD 

The kinetic equation for the normal Fermi liquid consisting of charged particles in the 
presence of compensating background was considered in [15] in the quasiparticle approxi- 
mation for the spectral function. It was shown that the spin splitting of the energy levels 
being neglected, the spectral function for a system in a nonquantizing magnetic field can be 
written in the form: 

r(pu; RT) 



(28) a(pu;RT) 



[u> - ^f RT » 2 - U {RT) - a"F(p; RT) - Re a(pu; RT)} 2 + 



2m 

where p is the canonical momentum and u(RT), A(RT) are scalar and vector potentials of the 
electromagnetic field, A(RT) being chosen in the Coloumb gauge divA(RT)=0. The kinetic 
equation for quasiparticle distribution is written in the form (25) [T5] and acquires a usual 
form after the transition to kinetic momentum P = p — (e/c)A(RT) [lj. We would like to 
stress that the kinetic equations for normal Fermi liquid were written in [TO] right in the form, 
corresponding the expression (28) for the spectral function, concerning the dependence of all 
the quantities on vector potential A and canonical momentum p and then were transformed 
into gauge-invariant form [T7]. More details about the development of transport theory of 
interacting fermions in an electromagnetic field can be found in [T8]-|20j. 

The derivation of the phenomenological Landau-Silin kinetic equations in the case when 
spin splitting of energy levels is taken into account was produced in [21] also in the quasipar- 
ticle approximation. In this case all quantities become matrices in spin space, in particular 
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the quantity e(pu; RT) has the form: 

r(ft — -A(HT)) 2 i A 

(29) e(pu; RT) = — — £— i — + u(RT) I S .curlA(RT) + a(pu; RT), 

l 2m J mc 

A 

where I and S are the unit matrix and the set of three spin matrices [1J, a is a Hermitian part 
of the self-energy matrix a. The spectral function is still given by Eq. (28). The spectral 
function a and the inverse life-time of a particle's state T are Hermitian matrices in the case 
under consideration. 

When T is considered to be finite (nonzero), the next possibilities can occur. If the 
magnetic field is not strong, so that T exceeds the spin splitting of the energy levels, the last 
is not essential and e(pu; RT) can be written in the form: 

(30) e(pu; RT) = ^~ c A ( RT )) + + .RT) + Re a(pu; RT). 

2m 

The approximation for the spectral function in this case can be written in a way analogous 
to Eq. (20): 

2 ~ r 



(31) a(pu; RT) = 2nZ(p; RT)5(lu - E(p; RT)) + Z 2 P 



(c^-S) 2 ' 

where E is a solution of Eq. (10) with e given by Eq. (30). 

The situation turns to be more complicated when spin splitting of the energy levels should 
be taken into account. More convenient form of the theory suitable for the generalization of 
the quasiparticle approximation for the spectral function, than that presented in [21], was 
developed in [22]. First of all , we should stress, that even in the quasiparticle approximation 
for the spectral function the kinetic equation of the phenomenological theory are valid only 
with the precision to the squared ratio of spin splitting of the energy levels to the chemical 
potential of the system [22]. Thus, producing the derivation of the kinetic equation with 
the improved approximation for the spectral function we should take into account only the 
terms that do not exceed this precision. 

We start with the expansion of the matrices g < and e over the full set of matrices in spin 
space, 

(32) g< = ig / + 5.g, 
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t A 

(33) e = eJ + 2S.e 2 . 

These expansions should be substituted into the generalized KB kinetic equations which 
under consideration are written in the form [22]: 

(34) (ul-e+ ^r)g< + b -[ul - e + ^T, g<] - a<(g + L -a) - ^[a<, g + L -a\ = 0, 

(35) g<(ul -e- L -T)+ L -[g<M-e- L -T]-{g- L -a)a< - L -[g - l -a, o<\ = 0, 
where g is a Hermitian part of the matrix g(pZ; RT). 

Separating in Eqs. (34) and (35) the terms corresponding to the collisionless case, we get in 
the quasiparticle approximation (r=0) the equations: 



(36) [u - ei,g ] - [e 2 i,gi] = 0, 

(37) [u- ei,|| - [e 2) g ] -2e 2 x |= 0, 

(38) (u - ei)g - e 2 -g= 0, 

(39) (w - ei)gj - e 2i g = 0. 



Here, the quantities with Latin indexes are the cartesian projections of the corresponding 
vectors; summation over the repeated indexes is supposed to be done. Eqs. (36)-(39) lead to 
the kinetic equations of the phenomenological Landau-Silin theory [22] in the limit r — > 0: 
Eqs. (38) and (39) give the expressions for the spectral functions, while the Eqs. (36) and 
(37) give kinetic equations on the basis of the determined spectral functions. 
With the help of Eqs. (38) and (39) we get the following expressions for the functions g 
and g in the quasiparticle approximations [22J: 



(40) g = / T 5(a;-e T ) + / i 5(a;-e i ), 

(41) g = fa(u - e T ) + faiu - e ; ) + fS(u - ei ), 
where 

(42) e T = ei - \e 2 \ ; e l = e 1 + \e 2 \. 



Vector fi is antiparallel to e 2) fi is parallel to e 2) and / is perpendicular to e 2 . 
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The argument of all the delta functions in (40), (41) except one of them, are not equal 
to u — e\. Thus, the substitution of expressions (40), (41) into the Eqs. (36) and (37) will 
generate the terms with the derivatives of the 5— functions. As a result, after rather tedious 
transformations of Eq. (36) we get the Eq. (19.41) from [22]. The analysis of the distinction 
between different renormalizing factors entering this equation, leads to the mentioned above 
conclusion about the precision of the validity of the phenomenological kinetic equation in 
the case T = 0, i.e at zero temperature. In the case of finite T, (but T being less then the 
spin splitting) it is necessary to put down the system of Eqs. (36)-(39) in the collisionless 
approximation, but saving the terms with T, which do not enter the collision integrals. It 
can be done, tracing carefully what terms with T in the equations in the absence of magnetic 
field do not enter collision integrals. Thus, we come to the system of equations: 

(43) {ul - e)g< + L -[uI- e, g<] - a<g - i[a<, g] = 0, 

(44) g<(cul - e) + L -[g<, ul - e] - ~ga< - ±[g, a<] = 0, 

where the expansions of the type (32), (33) should be done for the functions g and a < = Tf: 

1 A 

(45) g =-g~ / + S\f, 

A 

(46) a< = a 1 < I + 2S.a 2 < . 

Finally, taking into account the comment about the terms, exceeding the precision of the 
equation's validity, we get the system of equations that will lead to the kinetic equations: 

(47) [u - ex, g ] - [e 2i , gj + [g , <r<] = 0, 

(48) [u - ei , H - [e 2 , - 2e 2 x g + [g , <7 2 <] = 0, 

and the system of equations that will give the expressions for spectral functions: 

(49) (u - ei)g - e 2 .g - g Q a< = 0, 



(50) 



- e i)& - e 2ig - So^l = 0. 
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Now we consider the case e-p = ei — Je^l and choose the approximations for g and g in the 
form: 



(51) g = Z^5{u -E^) + Z 



2 °f 



(uj-E) 2 ' 



(52) g = Z t / t 5(cj - £ T ) + Z 2 ° 2 



where and E are the solutions of the equations: 

(53) E(p;RT) = e 1 (p;RT)l =E , 

(54) £ T (p; RT) = [e a (pk; RT) - \e 2 (p; RT)\}l =E , 
and renormalizing factors Z^ and Z are defined by the expressions: 

(55) Z i \p ] RT) = l- d - e ^* T) 



du 

d e^pu-RT) 



u)=Et 



(56) Z-\p-RT) = l 

If we substitute the expressions (51) and (52) into Eqs. (49) and (50), these equations would 
be satisfied with the precision to the terms of the order T 2 , |e~2| 2 and r | | - The substitution 
of expressions (51), (52) into Eqs. (47) and (48) lead to the Landau-Silin kinetic equations 
with the precision to the terms of the indicated order. The same is valid, if we consider the 
case ej_ = e± + | | - Thus, the kinetic equations of the phenomenological theory, (the Eqs. 
(7.21) and (7.23) in [22]) are valid in the case of finite temperature up to the terms linear in 
T, provided T does not exceed the spin splitting of the energy levels. 
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